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ABSTRACT 

Recent adaptive optics systems in astronomy achieve high-angular resolution. With 
the extreme stability of the images, detection at very low fluxes can be reached using a 
coronograph at the diffraction limit of the telescopes. This paper is an overview of the 
issues of stellar coronography used at the diffraction limit. Image formation through 
such a system is illustrated by numerical simulations. The description of a coronograph 
implemented on the VLT adaptive optics prototype, COME-ON, is presented as well 
as the first observations. 

Subject headings: instrumentation: miscellaneous - techniques: miscellaneous 
- stars: imaging - atmospheric effects 

1. Introduction 

Until the 1930's, the solar corona was only observed during total solar eclipses. To avoid this 
restriction, Lyot invented the solar coronograph (1930, 1931, 1939). This instrument was a solar 
telescope equipped with a mask, whose size was exactly that of the sun. The light coming from the 
sun disk was blocked, which allowed the corona, one million times fainter than the photosphere, 
to shine through. Later, the light rejection rate was improved by modifying slightly the original 
layout (Evans 1948; Newkirk & Bohlin 1963; Fort et al. 1978; Koutchmy & Belmahdi 1987). 

Recently astronomers have tried to apply this method to observe stellar environments. By 
occulting the stellar luminous contribution, the circumstellar material becomes detectable. In 1984 
a coronograph made of a reticle in the focal plane enabled astronomers to detect a circumstellar 
disk around the star f3 Pictoris (Smith & Terrile 1984, Vilas & Smith, 1987). More recently, an 
infrared coronograph has been designed in order to detect brown dwarfs (Beuzit et al. 1991). 
The limiting factors for this kind of stellar coronograph are the lack of spatial resolution of the 
instruments and the very faint brightness of the stellar surroundings. To observe a planet separated 
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by less than 5 AU from a star located at 5 pc, spatial resolution better than 1" is required. This is 
also the resolving power limit of present telescopes. To achieve a good efficiency, precise centering 
of the star on the coronographic mask is necessary. Due to atmospheric disturbance, this cannot 
be done efficiently unless the mask is large. Thus Clampin et al. (1992) have built an optimized 
coronograph suited to the good imaging of the ESO NTT telescope at La Silla (Chile). 

A clear way to avoid atmospheric turbulence is to perform space-based observations. Following 
this idea, Bonneau et al. (1975) proposed a coronographic mode for the Faint Object Camera 
(FOC) of the Hubble Space Telescope. Breckinridge et al. (1982) also proposed a coronographic 
camera for the Space Telescope and Watson et al. (1991) suggested the use of a coronograph on a 
large aperture lunar-based telescope. 

Another way to improve image quality is to use active optics techniques. The "Shift & Add" 
method used in real-time with a tip-tilt mirror, as developed by Durrance & Clampin (1989) 
at John Hopkins University then by Golimowski et al. (1992) and also for the camera HRCam 
at the Canada- France-Hawaii telescope (Racine & McClure 1989), improves the efficiency of a 
coronograph and enables reduction of the occulting mask diameter. Current adaptive optics 
systems can also be used. This is the approach chosen that will be described hereafter in this 
paper. With recent adaptive optics systems (Rigaut et al. 1991; Beckers 1993), the resolving power 
is limited by the telescope diffraction, less than 0.2" in the near-infrared, and the resulting images 
are motionless. The optical star beam can be centered on the occulting mask with high precision 
and is available for standard coronographic processing. This design has been implemented in 
January and April 1991 on COME-ON (Malbet 1992a, 1992b), the VLT adaptive optics prototype 
at the ESO 3.6 meter telescope in La Silla (Chile). This system and the first observations are 
described in this paper. The paper also addresses the performance of the first version of this 
coronograph. In the summer of 1991, the prototype COME-ON was updated to COME-ON+: 
the optical scheme was changed and the coronograph removed. In parallel, Durrance & Clampin 
(1989) have designed a full adaptive optics system dedicated to coronographic imaging at John 
Hopkins University. Unfortunately the full system has not been built to date. 

In order to estimate the science that can be achieved with high angular resolution 
coronography, the reader is referred to some representative papers published in classical stellar 
coronography: imaging of f3 Pictoris disk (Smith & Terrile 1984; Paresce & Burrows 1987), 
circumstellar imaging of symbiotic stars, planetary nebulae (Paresce, Burrows & Home 1988; 
Burgarella, Clampin & Paresce 1991; Burgarella & Paresce 1991, Robberto et al. 1993) and 
luminous blue variables (Paresce & Nota 1989; Burgarella & Paresce 1991), imaging of the 
supernova SN 1987A (Paresce & Burrows 1989). With the JHU Adaptative Optics Coronograph, 
stellar coronography has improved in spatial resolution, e.g. for imaging blue luminous variables 
(Nota et al. 1992; Clampin et al. 1993). Images with very high contrast have also been obtained 
around young stars (Nakajima & Golimowski 1995) and a very-low mass star has been detected 
around a binary (Golimowski et al. 1995). These results show the potential of a coronograph 
used on a full adaptive system. Recent results with the new coronograph of COME-ON+ by 
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Fig. 1. — General optical schematic of a Lyot coronograph. Four optical planes are represented. 
Optical rays have been traced from 5 different points on the entrance pupil plane and for 3 different 
field rays. 

Beuzit, Mouillet & Lagrange (1995) include imaging f3 Pictoris disk up to 1" from the star in the 
near-infrared. 

Image formation through a coronograph is treated in Sect. 2. Numerical simulations illustrate 
the performances in Sect. 3. Description of the coronographic mode of COME-ON and the first 
results constitute Sect. 4. 



2. Image formation through a coronograph 



The simplest coronograph is composed of an occulting mask introduced in the optical beam 
at the focus of the telescope. Lyot (1930, 1931) pointed out that the efficiency of a coronograph 
is appreciably increased, after the focal mask, by the use of a diaphragm in a pupil plane slightly 
smaller than the telescope aperture. A coronograph is then made of a focal occulting mask and a 
pupil stop, usually called Lyot stop. The mask is optically conjugated with the detector, and the 
Lyot stop with the telescope pupil (see Fig. 1). 

In this section, I assume the adaptive optics system is perfect and restores a perfect wavefront. 
Therefore, I will consider the image formation in the far field approximation. Bonneau et al. (1975) 
and Labeyrie (1985) provided good presentations of this process, but without any mathematical 
formulation. Noll (1973) quantified the contribution of the Lyot stop for the intensity calculated 
on the optical axis. Wang & Vaughan (1988) have proposed an approximation to obtain the 
intensity in the focal plane in cases where a Lyot stop is used. The aim of this section is to describe 
the image formation through a coronographic system within the Fraunhofer diffraction theory, 
where the lenses are assumed perfect. In this section, I will show that the coronographic process 
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is similar to a convolution, but with a point spread function different for each point of the sky. 

A coronograph consists of cascaded optics. There is first an occulting mask which has two 
purposes: 

• to reduce the light coming from the central star by a direct occultation 

• to filter out the light at low spatial frequency. 

The results are twofold: 

• the on-axis stellar light in the following pupil plane is almost totally cancelled, depending on 
the size of the occulting mask, 

• some light passes through. As the mask removes the low spatial frequencies in the image 
plane, the remaining light is on the edge of the pupil which corresponds to the highest 
frequencies. 

The Lyot stop blocks the light located at the pupil edge. For a point located far from the central 
star, the occulting mask has a weak effect and the Lyot stop decreases the size of the aperture 
only slightly. Since the Lyot stop is located at a pupil plane, the removal of stellar light takes 
place everywhere in the focal plane, not only in the central part as with the occulting mask. In 
a way, the occulting mask is useful to avoid saturation of the detector with central stellar light, 
but not useful for detecting circumstellar material, because the diffraction pattern remains at the 
same level of intensity outside of the occulting mask. The purpose of the Lyot stop is to decrease 
this level of diffraction. 

2.1. Instantaneous amplitudes 

The amplitude of light arriving onto the focal mask is the convolution of the initial light 
amplitude of the object, tp (Q), with the normalized impulse response of the telescope, Gy(0): 

= g t (0)*</> o (0), (i) 

with = (6 X , Oy) is the position vector in terms of angular size. 

The field amplitude just after the occulting mask is -i/> m (0) multiplied by 11^(0), the 
transmission function of the mask. Therefore the resulting amplitude on the detector is: 

</>(©) =G L (0)*{lI M (0) [g t (©)*</> (©)]}, (2) 

where G^(0) is the impulse response associated with the Lyot stop. By permuting the integrals 
of the two convolution operations, one finally gets: 

</,(©) = J rf0 ^/>o(0 )<?(©; 0), (3) 
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The amplitude in the following pupil plane just before the Lyot stop, is the Fourier transform 
of the amplitude found just after the occulting mask: 



^(/) = UM(/)*[n T (/)^ (/)], 



(5) 



where LI^ is the pupil function (the Fourier transform of the impulse response Gt), Um is the 
Fourier transform of the mask and / = (u, v) is the spatial frequency coordinate. It can be written 
as: 

W) = / rfeVo(e)/c(e;/), (6) 

where 

K{Q-J) = J d@U M (@)G T (@-@)e-' 2 ^ f. (7) 



2.2. Intensity distribution 



The intensity is defined as / = (tptp*). In the case where the incident light is coherent, 

(M®iW &)) = M®iW &), (8) 

one gets respectively in the focal plane and in the Lyot plane: 

/(e) = ^(e)^(e), (9) 

with V(6) = J deVo(© )£(©;©), (10) 

i L (f) = ^(/> L m (ii) 

with i, L {f ) = [ deVo(e')/c(e ;/). (12) 



(13) 



In the case where the incident light is incoherent, 

(^o(e 1 )^(e 2 )) = /o(e 1 )5(e 1 -e 2 ), 



where Io(@) is the intensity distribution of the object and delta(&) is the delta-function. One gets 
respectively in the focal plane and in the Lyot plane: 



/(e ; 



i L (f] 



dO 7 (e )J r (@;@) 
with T(Q-Q ) = Q(Q ■& ) 

dQI (Q)U(Q-J) 
with %(&■ f)= JC(@-J ' 2 



(14) 
(15) 
(16) 
(17) 
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2.3. Point Spread Function of a coronograph 

The notion of point-spread-function (PSF) used in traditional optics is not valid in the context 
of a coronograph. The resulting intensity is not the convolution of the incident intensity by a PSF, 
because the modulation transfer function is no longer invariant to translation. An unresolved 
star will have a different response when located under the occulting mask than that it has when 
located outside the mask. However that is the intent of a coronographic system, to enhance low 
flux level while attenuating the unwanted axial light. 

. . ->' 

If the object is an unresolved star at position , the intensity distribution is 
/ o (0) = I S(Q - ), then 

= I T{Q ;0) (18) 

Therefore J-{® ; ) is the "point spread function" associated with the direction . There is not 
a unique PSF but a PSF for every point on the object. This set of PSFs, J-(@ ;0), is easily 
calculated if one notices that it is the image of a point source at direction through the system. 
Therefore J-(@ ; ) is only the PSF of the telescope occulted by the focal mask and convolved 
with the PSF of the Lyot stop. 

In terms of discrete sampling, if the position on the detector is identified as one gets: 

= ~^jIok,iFk,i,ij, (19) 

k,l 

where F^ n j is the 4D-tensor associated with the system transfer function. To recover the original 
7 , the tensor F has to be known and inverted. 

The shape of 7i(@ ;0) can be understood by the effect of diffraction of a star centered on 
the focal mask. This effect is to convolve the light distribution in the pupil plane by the mask 
impulse response (Bonneau et al. 1975; Labeyrie 1985). Since only the high spatial frequencies are 
not filtered by the occulting mask, the light is mainly suppressed far from the pupil edge. This 
results in a pupil which is over-intensified near the edges. For a star not centered on the optical 
axis, the focal mask has almost no effect. Therefore the light distribution in the pupil plane is not 
changed. The Lyot stop suppresses the bright edge fringe from the central bright star, but has 
minor consequences for an object far from the mask. The Lyot stop must be sized to the mask 
diameter so that the blocked area has a width equal to that of the mask impulse response in the 
pupil Lyot plane. 

The coronograph PSF beyond the focal mask edge is very similar to the Lyot stop PSF. The 
Lyot stop PSF is obtained if the entrance pupil were the Lyot stop. Likewise the coronograph PSF 
far from the mask is obtained by removing the focal mask and not the Lyot Stop. In a complete 
coronograph, the detector intensity distribution is the convolution of (1) the intensity distribution 
just before the mask multiplied by the mask , and (2) the Lyot Stop PSF. When the off-axis 
direction is far from the mask edge by more than the width of the Lyot stop PSF, the intensity 
distribution is very close to the one obtained with no focal mask, but with a pupil. 
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Fig. 2. — Intensity distribution in a coronographic system for detection of a binary star with a 
flux ratio of 1000 and a separation corresponding to the radius of about 20 bright Airy rings from 
the center. The focal mask stops the light inside 12 Airy rings. The panels show the intensity 
distributions: (a) in the entrance pupil plane, (b) in the first focal plane before the occulting mask, 
(c) in the first focal plane after the coronographic mask, (d) before the Lyot stop, (e) after the Lyot 
stop and (f) in the final focal plane of the detector. The logarithmic color scale is identical for the 
different images of the focal planes and also for the images of the pupil planes. 

Figure 2 shows the different intensity distributions along the optical path of a coronograph in 
the case of a close binary. Note the over-intensity at the pupil edge in Panel (d) which has been 
removed in Panel (e) by using of a Lyot stop. 



2.4. Using the Lyot stop 



The Lyot stop decreases the intensity of the stellar diffraction wings at every point in the 
focal plane even far from the occulting mask. It modifies only slightly the intensity of an object 
located outside the coronographic mask. Therefore the contrast between the off-axis object and 
the level of intensity from the diffraction wings of the on-axis star increases. Without a Lyot stop 
the level of intensity of the stellar diffraction wings outside the occulting mask is not changed 



- 8 - 




Fig. 3. — Intensity in the Lyot pupil plane for two kinds of coronograph: (a) occulting mask with 
a same size as the resolving power (1 Airy disk); (b) occulting mask size large compared to the 
resolving power (12 Airy rings). The figures display the intensity distribution in the pupil plane 
following the occulting mask. Note that the over-intensity due to the diffraction of the light on the 
edge of the occulting mask is sharper in panel (b) than in panel (a). 

by the coronograph and the contrast between the object and the background level remains the 
same as it would be with no coronographic device. The contrast enhancement comes from the 
attenuation of the edge discontinuity in the pupil plane by the Lyot stop. This discontinuity size 
is inversely proportional to the occulting mask size (Fourier transform properties). Figure 3 shows 
the difference between a coronograph with a mask size large compared to the telescope PSF size 
(panel a) vs. a coronograph whose mask size is similar to the PSF size (panel b). 

• If the mask size is large compared to the full width at half maximum (FWHM) of the 
telescope impulse response (Panel (a) of Fig. 3), then the over-intensity in the Lyot pupil 
plane will be sharp and may be suppressed by a Lyot stop slightly smaller than the pupil. 
This is the case for solar coronographs and also for stellar coronographs that use large masks 
compared to the FHWM of the image. 

• If the mask size is of the same order of magnitude as the FWHM of the telescope PSF (Panel 
(b) of Fig. 3), then the over-intensity width is of the same order of magnitude as the pupil 
diameter. It is not possible to suppress it without significantly decreasing the Lyot stop size, 
resulting in a loss in resolution and collecting area. This is the case for coronographs used 
at high angular resolution (i.e. at the diffraction limit of the telescope), and the Lyot stop is 
not recommended. 
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If one is interested in the close neighborhood of an object, a small mask must be used. The 
Lyot stop will not significantly decrease the intensity of the central object. However if one is 
interested in objects further from the central object, the larger the focal mask, the sharper the 
over-intensity in the Lyot plane and the better the rejection rate. 



2.5. Summary 

Image formation occurs similarly in a coronograph as in a simple optical instrument. However, 
the coronographic mask breaks up the system transfer function invariance. Therefore the point 
spread function depends on the object position in the field. The convolution operation is then a 
multiplication by a 4-fold tensor - the transfer function of the coronographic system. However, at 
a distance corresponding to the full width at half maximum (FWHM) of the Lyot stop impulse 
function from the focal mask, the coronograph PSF is similar to the Lyot stop PSF. 

If the size of the coronographic mask is large compared to the telescope PSF, the coronograph 
has outstanding performance. In that case, over-intensity at the Lyot pupil edge is sharp and is 
easily stopped by a Lyot stop few percent smaller than the telescope entrance pupil. Obviously a 
coronograph works better with large masks. On the other hand, stellar coronography and studies 
of the circumstellar environment demand observations close to the star. Therefore very high 
angular resolution combined with coronography is required, but one must use a smaller mask 
(whose size is comparable to the FWHM of the PSF) and the Lyot stop is not as efficient. 

The following section presents numerical simulations. Coronographic performances will be 
discussed in terms of the mask and Lyot stop size versus cutoff frequency of the telescope. 



3. Expected performances 

In this section, I apply previous results to a simple case, and compute numerically the 
intensity at the Lyot plane, I L (f), and at the focal plane /(©). This section demonstrates the 
dependence of these quantities on both mask and Lyot stop size. 

Consider the case of a circular pupil with cutoff frequency f c , a circular occulting mask 
with angular diameter 6 m , and a circular Lyot stop with cutoff frequency af c (with a > 0). 
Furthermore, suppose the magnification between successive planes is one. The transfer functions 
and impulse responses used in the previous section are found in Table 1. As described at the 
end of Sect. 2.4., a should be such that the over-intensity in the Lyot pupil plane is suppressed. 
The diameter of the dark ring closest from the pupil edge has been chosen for the Lyot Stop 
diameter. Therefore there is a relation between computed for occulting mask diameter 6* ml and 
a> 2 computed for occulting mask diameter m2 : 

a2 = l_(l_ ai )^l. (20) 
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Table 1: Mathematical expressions for the transfer functions and impulse responses of the telescope 
entrance pupil, the mask plane, and the Lyot plane in the case of a circular pupil. 



Location 


Direct Planes (focus) 


Fourier 


Planes (pupil) 


Telescope aperture 


G T (e) = Ji(>r/.«)/v/s> 


n T (/) : 


f 2/v^/ c if p < f c /2 
{ if p > f e /2 


Mask 


u M( e /(U = (» 

{ 1 ' \ 1 if 6 > e m 




= s(f)-e m J 1 (7r0 m p)/2p 


Lyot stop function 


G L @) = af c J 1 (waf c 0)/20 




f 1 if p < af c /2 
\ if p > af c /2 



with 0= ||©|| = + and p = Vu 2 + v 2 



The optical analysis software, COMP (Controlled Optics Modeling Package, Redding et al. 1992), 
was used to calculate the intensity distribution at different planes in the coronograph. 

One may wonder whether results computed with simple assumptions (circular pupil and 
circular Lyot stop) are useful for understanding actual systems, since telescopes usually have 
secondary obscuration and spider? The final performance is slightly different, but the effects are 
basically identical: the main effect comes from the diffraction of light at the edges of the pupil, 
either at the outer edge of the pupil, at the edge of the secondary obscuration or at the edge of 
the spider. A circular pupil makes computing but also interpretation of the results easier. 

3.1. Impulse responses of a coronograph at the Lyot stop and focal planes 

Figure 4 and Fig. 5 display the functions 7i(@ ; f ) and J-(® ;©). A squared mesh of 
1024 X 1024 pixels was used. The FWHM is equal to 32 pixels (l//c). The PSFs are computed for 
two mask diameters, 39 pixels corresponding to the first dark ring and 104 pixels corresponding 
to the third dark ring of the telescope PSF. The radial variations of J-(@ ;©) and 7i(@ ;f ) are 
plotted for three PSF positions corresponding to: 

• the optical axis (Panels a and d) 

• the third bright ring of the Airy pattern centered on the optical axis (Panels b and e) 

• the fifth bright ring of the Airy pattern centered on the optical axis (Panels c and /) 

The effect of the coronograph is demonstrated in Panels a and d of Fig. 4 and Fig. 5. In Fig. 4, 
one can see the over-intensity surrounding a dark pupil. It is not the case for the off-axis directions 
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Fig. 4. — Pupil plane impulse response H(Q ;/) of a coronograph with a circular entrance pupil. 
The scale is identical in the 6 panels. 

(a), (b), (c) corresponds to a coronograph whose focal mask has the same diameter as the first Airy 
dark ring, (c), (d), (e) belongs to a coronograph whose focal mask has the same diameter as the 
third dark Airy ring. 

(a) and (d) are the intensity distributions in the Lyot plane for 0=0. (b) and (e) are the intensity 
distributions in the Lyot plane for corresponding to the third dark Airy ring, (c) and (e) are 
the intensity distributions in the Lyot plane for corresponding to the fifth dark Airy ring. 



(Panels b, c, e, /). The difference between Panel a and d is the width of the over-intensity due to 
the size of the occulting mask. In Fig. 5, the influence of the size of the occulting mask. First the 
spatial resolution is less in upper panels compared to lower panels. Secondly the on-axis direction 
does not present a volcano shape and the central peak is dominant. However with both occulting 
mask sizes, the gain is evident in intensity. The color scale is logarithmic and the intensity for 
the off- axis direction is much more bright than for the on-axis direction. One can also notice that 
the effect of the mask is almost canceled for the directions far from the mask edge (Panels b, c, /). 
Quantitative results are given in the following section. 



3.2. Rejection rate and intensity gain 
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Fig. 5. — Detector plane impulse response J-(Q ; / ) of a coronograph with a Lyot stop (a ^ oo). 
The scale is identical in the 6 panels. 

(a), (b), (c) corresponds to a coronograph whose focal mask has the same diameter as the first Airy 
dark ring, (c), (d), (e) belongs to a coronograph whose focal mask has the same diameter as the 
third dark Airy ring. 

(a) and (d) are the intensity distributions in the focal plane for 0=0. (b) and (e) are the intensity 
distributions in the focal plane for corresponding to the third dark Airy ring, (c) and (e) are 
the intensity distributions in the focal plane for corresponding to the fifth dark Airy ring. 

The rejection rate R of a coronograph is by definition 

R=?j^, (21) 

where I w is the total intensity of the star in presence of coronograph and I w / is the total intensity 
of the star in absence of coronograph. It depends on the mask and Lyot stop sizes as well as 
the angular separation between the central star and the interesting object. Figure 6 displays the 
rejection rate as a function of the distance from the optical axis for 3 different mask diameters 
(first, third and fifth dark rings). The mesh is the same as the one used in previous section. 

Table 2 gives the result of the asymptotic rejection rate for a coronograph equipped with Lyot 
stop and a coronograph without one. Improvement is clear between the two designs. However the 
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Table 2: Expected rejection rates of a stellar coronograph at high spatial resolution 



Mask diameter 


with Lyot stop 


w/o Lyot stop 


Am 


1 rinp' 


20 


6 


2.0 


2 nd ring 


66 


11 


2.6 


3 rd ring 


120 


16 


2.9 


4 tn ring 


180 


21 


3.1 


th 

5 ring 


250 


26 


3.2 


10 tn ring 


580 


52 


3.4 



gain is less when the mask is small because the Lyot stop suppresses a large part of the pupil. 
The third column gives an estimate of the gain in magnitude on the stellar profile due to the 
addition of the Lyot stop. Am is equal to 2.5 \og(2R w / R w / ) and is the same everywhere in the 
focal plane, except near the occulting mask. The factor 2 comes from the new repartition of the 
light after the convolution due to the Lyot mask. About 50% of the light is diffracted in the rings 
inside the mask area. The numbers can be compared with the actual results obtained on the sky 
and presented in the next section. 

3.3. Expected limitations from imperfect actual systems 

Adaptive optics systems are experimentally not perfect. They only correct a certain number 
of modes. For example, tip-tilt compensation systems deal with only two modes. The COME-ON 
system, described in the next section, corrects 21 modes and the VLT adaptive optics is planned 
to correct about 256 modes. Not all the modes are useful depending on the wavelength and the 
telescope diameter. At A = 10/im, only tip-tilt correction is necessary for a 4-m telescope in a 
good site. With 21 actuators on a 3.6-m telescope, COME-ON provides almost full correction for 
wavelengths larger than the near-infrared K band (Rigaut et al. 1991). K band images themselves 
are not fully diffraction-limited. The limit between full correction and partial correction is given 
by Lord Rayleigh's criterion (Marechal & Francon 1970; Mariotti 1988): residual phase errors 
must be smaller than A/4. 

In the case of full correction, aberrations are small phase errors (A<^> <C 7r / 2) and can be 
understood as scatter-type aberrations. The phase contribution e 8A< ^ becomes 1 + iA(f> and the 
scattered light is linearly added to the image. The intensity of the scattered light is proportional 
to the source intensity in the image. The pattern corresponds to the Fourier transform of the 
phase errors convolved with the PSF. This phenomenon is similar to the scattered light which 
comes from the collecting mirrors of large telescopes and which limits the dynamical range of 
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1st 2nd 3rd 4th 5th 6th 7th 8th 9th 

Distance from the center (dark Airy rings) 



Fig. 6. — Coronograph rejection rate with Lyot stop. The diameter of the different masks are 
respectively: first dark ring (solid line), third dark ring (dashed line) and fifth dark ring (dotted 
line). The a;-axis represents directions at different angular distance from the optical axis. The size 
of the Lyot stop has been computed to match the first dark ring at the pupil edge. 

spatial telescopes (Brown & Burrow 1990; Malbet, Yu & Shao 1995). The coronographic image is 
then photon-noise limited by the scattered light and many exposures are needed to decrease the 
background noise by subtraction. This limitation does not originate in the coronograph, but in 
the adaptive optics system correction. 

At shorter wavelengths, adaptive optics systems work in partial correction and are no longer 
diffraction-limited. Image formation becomes much more complex. We try to provide some 
answers but the understanding of this phenomenon is still limited. Beckers & Goad (1987) have 
shown that long-exposure stellar profiles in partial correction are spiked helmet-shaped in which 
an Airy shaped spike is superposed on a broad halo-like background with a width approximately 
equal to that of the original seeing disk. This was confirmed experimentally by Rigaut et al. 
(1992). As a matter of fact, what happens in the image formation through a coronograph is 
difficult to understand without extensive simulations and experimental tests. Qualitatively, one 
can estimate that the Airy shaped spike behaves like the fully corrected image and the broad 
halo like the seeing-limited image. The main limitation would come from the halo as in classical 
coronography. 
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Table 3: Rejection rate of coronographs with no accurate mask centering. 



Mask diameter 
(centering error) 


with Lyot stop 


without Lyot stop 


{Om/W) 


(0 m /2) 


(0 m /lO) 


(0 m /2) 


l s ^ ring 


17 


15 


6 


3 


2 nd ring 


66 


42 


11 


7 


3 rd ring 


120 


47 


16 


9 


th 

5 ring 


250 


49 


21 


12 


10 tn ring 


580 


54 


52 


17 



An interesting case is when the star is not accurately centered on the coronographic mask. 
This is relatively usual since in actual systems the star is never very accurately centered on the 
mask. Let us assume a cr-wide gaussian distribution of the centering probability and that the 
fluctuations of the rejection rate comes only from the centering of the central bright star. Table 3 
presents the numerical results of this model by using the same procedure as in Sect. 3.2. to 
calculate the rejection rate. Rejection rates from a bad centering of one tenth of the mask size 
are not too degraded. For a bad centering of one half mask size, the rejection rates decrease by a 
factor of 2 (without a Lyot stop) up to a factor of 10 (with a Lyot stop). This is the reason why 
tip-tilt compensation system gives better results than classical stellar coronography. However the 
performance of tip-tilt compensated systems is still limited by the width of the stellar profile. 

Another potential limitation is the isoplanatic angle. The isoplanatic angle is a radius of a 
circle on the sky over which atmospheric wavefront disturbances are practically identical. If the 
averaged altitude of the seeing layer is 6 km, the isoplanatic angle is about 60" at A = 2.2/im and 
about 3" at A = 0.5/im. High angular resolution coronography with diameter of telescope larger 
than 2 m have always a spatial resolution much smaller than the isoplanatic angle and is likely to 
be limited by this effect. 

4. The coronographic mode of COME-ON 
4.1. Description 

COME-ON is a prototype astronomical adaptive optics system for the VLT. It has been 
designed and built in France (Observatoire de Paris, ONERA, and Laserdot) in collaboration with 
the European Southern Observatory (ESO). The instrument design and performance have been 
described by Rousset et al. (1991) and Rigaut et al. (1991). 

The COME-ON optical design permits the insertion of a stellar coronograph. The stellar 
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Fig. 7. — Optical layout of COME-ON coronographic mode. The solid lines represent the envelope 
of the optical beam when the coronograph is inserted. The dashed lines represent the original 
optical path of the light coming from COME-ON when the coronograph is not used. 

beam is split by a dichroic into a visible beam used for wavefront sensing and an infrared beam 
used for science. This infrared beam is corrected by the deformable mirror located before the 
dichroic beam-splitter. The outcoming beam is corrected for the atmospheric turbulence and the 
images are very stable. The coronograph design is inserted between the last pupil plane and the 
focus. This design allows an easy implementation for astronomical use. 

4-1.1. Optical layout 

The optical layout is shown in Fig. 7, where only the imaging arm of COME-ON is drawn. 
In order to implement all coronographic features, focal and pupil planes are created for the focal 
mask OM and the pupil Lyot stop LS. Three lenses LI, L2 and L3 are used to create these planes: 

• LI reduces the focal distance and forms the star image on the occulting mask OM. The 
optical aperture is then f/20 instead of the original f/40. 

• OM is a glass plate that supports the occulting masks, 

• L2 creates a plane conjugate to the entrance pupil plane, 

• LS is a glass plate that supports the Lyot stop, 
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• L3 forms a convergent beam identical to the original beam and images the masked star onto 
the detector. The optical aperture is then f/40 or f/80 depending on the camera optical 
configuration. 

These 5 optical elements are installed on an optical plate inserted between the On-Off mirror 
of COME-ON (M6 in the quoted publications) and the infrared detector. The M7 plane mirror 
only changes the beam direction to optimize space utilization. The lenses and glass plates which 
support the mask and Lyot stop are made from insoluble Fluoride glass (AFG 450 made by Le 
Verre Fluore) which transmits infrared wavelengths. 

4-1-2. Occulting masks and Lyot stop 

The masks were made using microphotolithographic techniques at LARCA, a laboratory of 
the Observatoire de Paris. The 5 step procedure is: 

• evaporation of chromium and then gold on the glass, 

• deposition of a photosensitive resin layer 

• ultraviolet illumination of the resin through a mask with the desired pattern, 

• chemical removal of the chromium and the gold in the open zones, 

• ultraviolet illumination of the rest of the glass support to remove the remaining resin. 

To allow for different types of observation, 8 masks of differing diameters were placed on 
the glass support. A remote controlled translation stage shifted from one mask to the next and 
positioned the mask over the star. The mask sizes were chosen to approximately fit the diameter 
of the first dark diffraction rings at the different infrared wavelengths J, H, K, L and M. The 
chosen sizes were: 

95, 120, 175, 205, 215, 250, 370, 460 fim 

corresponding on the sky to 0'.'27— 1'/3. 

The Lyot Stop was made with the same technique. Only one stop was used. It was designed 
proportionally to the pupil shape: a circular aperture with a secondary obscuration and 4 spider 
arms. The pupil was reduced by 10% in size in order to increase the rejection rate. The glass 
support was mounted on a translation stage allowing fine positioning of the stop. 
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Table 4: Elementary integration for the observed star and performance of the coronograph 



Elementary integration time 


Gain 


Gain 


Rejection 


with coron. w/o coron. 


in time 


in intensity 


rate 


Sirius 192 ms 16 ms 


12 


9.73 


51 



4.2. Observations 

Observations using the coronographic mode on COME-ON were conducted January 1991 
and April 1991 with the ESO 3.6-m telescope during nights dedicated to young stellar object 
observations. Only preliminary tests were performed to qualify the system for future use. 

Images were recorded with a 32 X 32 InSb IR camera array (Lacombe et al. 1989). The pixel 
scale was 0'.'056 on the sky using the f/80 camera configuration. Chopping allowed for alternating 
on-source and sky exposures. Standard dead pixel elimination, flat fielding, and specific cross-talk 
correction between lines were performed for each object and sky elementary frame. The sky was 
subtracted from each object elementary frame, and the frame's total intensity was normalized. 

One star was observed: Sirius. The images were taken in the infrared K band. This star 
is bright and the integration time was set at the detector saturation limit. Twenty elementary 
images were added together to create an image. The optical layout for the coronographic image 
used an occulting mask and a Lyot Stop. For the image without the coronograph, the occulting 
mask was taken out, but the Lyot stop remained. Elementary integration times with and without 
the coronograph are summarized in Table 4. The spatial achieved resolution has been estimated 
at 0.3". This resolution is not diffraction-limited on this system because at this wavelength 
the wavefront correction is only partial, but the resolution is much better than any previous 
coronographic observation. Quantitatively, Rigaut et al. (1991) gives a residual phase error of 
11 rd 2 for A = 0.5/im, which corresponds to an absolute phase error of 0.3/im or A/3.7 in the K 
band., slightly smaller than the Rayleigh's criterion. 

The size of the mask used was 175 /im corresponding to 0.5". It was a little smaller than the 
second dark Airy ring, and with the achieved resolution, occulted the central core of the PSF. The 
smallest coronographic mask used previously in other coronographs is 4" (Golimowski et al. 1993). 

4.3. Results 

In this section, I describe the specifics of coronographic reduction process and I provide the 
technical results of the observations. In addition of the calculation of the rejection rate, I show 
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Fig. 8. — Interpolated images of Sirius (a) without and (b) with the COME-ON coronograph. The 
elementary integration time has been multiplied by 12 between the two images to reach the same 
level of flux. The intensities have been normalized in function of the intensities in image(a). Image 
(c) is the image of the residuals from the subtraction of image (a) from image (b). The histogram 
shows that the residuals in image (c) are smaller than the background noise level in images (a) and 
(b). Sirius was close to the edge of the detector in image (a) producing the "white" vertical band 
in subtracted image (c). See text for detailed comments on data processing. 



that the coronograph does not destroy the star radial profiles beyond a radius equal to the sum of 
the mask radius and the FWHM. 

Upper panels of Fig. 8 show the images of Sirius taken with and without the coronograph. 
The dynamic range was about the same in both series of images but the elementary time was 
multiplied by a factor 12 for the coronographic data. This gives an estimate of 12 for the gain in 
intensity. A better estimate is given further. 

The coronographic reduction process consists in subtraction of the PSF from the coronographic 
image. In this paper the PSF is the same star observed with the coronograph but without the 
occulting mask (however the Lyot stop remains). Therefore I show that after accurate centering, 
subtraction of the PSF image from the coronographic image gives only background noise. The 
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different steps of the process are: 

• Linear interpolation of images (a) and (b) by a factor of 4 in order to get enough precision 
in the centering. 

• Computation of software mask by convolving a hard-edge circular mask of 0.5" centered on 
the centroid of image (a) with the PSF (image (a)). 

• Multiplication of the convolved mask by the hard-edge mask in order to suppress the weight 
of the central pixels in the cross-correlation. 

• Cross-correlation of image (a) multiplied by the software mask and image (b) in order to 
know the angular separation between images (a) and (b). This method has the advantage to 
use only pertinent information (the intensity distribution outside of the occulting mask) and 
is unaltered from the remaining light at the mask. 

• Translation of image (a) by a distance corresponding to the calculated separation. 

• Computation of the ratio a between the two images by a least-square method. The \ 2 
metric is given by: 



where a,j and bij are the elements of image (a) and image (b) padded with zeros and M,j are 
the elements of the software mask used for the cross-correlation and centered on the star. 

• Image (c) is the result of image (b) minus a times image (a), then multiplied by the software 
mask. 

The results obtained for Sirius are shown in Fig. 8. The coefficient a is 9.73 corresponding to 
the gain in intensity. The rejection rate for the star is calculated by dividing the flux in image (a) 
by the flux in image (b) and is 51. The rejection rate for Sirius is very close to the theoretical 
value calculated in Table 2 where the rejection rate for an obscuration two Airy rings wide is 
66 for a coronograph with a Lyot stop. The difference can be attributed to misalignment errors 
and to optical bench flexures. To show the performances of the system, a histogram of the three 
images has been computed in Fig. 8 (Panel d). One can see that the residuals (image (c)) are 
below the background noise floor of image (a) and image (b). This means that the residuals from 
the subtraction are at the noise floor. This noise floor is reached mainly because of the noise level 
in image (a). 

The reduction process for coronography at the diffraction limit described in this section can 
be adjusted, but one still needs calibration images of the PSF (used to smooth the hard-edge 
mask). If a better signal-to-noise ratio is desired, a coronographic image of the PSF is required in 
order to get a smaller noise level. Then the coronographic image of the PSF can be subtracted 
from the coronographic image of the object. 
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In order to convert this result in useful astronomical units one can say the result corresponds 
to a gain of 2.47 ± 0.02 in magnitude for a mask size of 0.5". At a distance of 1" from the center 
of the mask, the level of star diffraction is then 12.5 magnitude fainter than the peak of flux from 
the star. 

5. Conclusion 

This paper presents an overview of the coronographic process for high angular resolution 
imaging. It shows high angular resolution coronography is a powerful technique that can reveal 
circumstellar environments. 

The occulting mask and Lyot stop must be carefully chosen. A relationship exists between 
them and the instrument resolution. The larger the occulting mask is relative to the PSF, the 
smaller the pupil reduction of the Lyot stop. If the mask is close to the PSF size, a Lyot stop 
is not necessarily required for it does not improve the rejection rate and only reduces spatial 
resolution. Simulations show the coronographic impulse functions in different planes and illustrate 
the coronograph principle. In the case of bad centering, the image data are still valuable, but the 
dynamic range and the rejection rate are lower. Therefore accurate centering is important for 
efficient observations and a better signal-to-noise ratio, but does not degrade data. A coronograph 
mode has been implemented on the adaptive optics prototype COME-ON. Initial observations 
indicate the coronograph has the expected behavior and illustrate the potential of such a system 
for astronomical observations. 

I would like to thank the COME-ON team, the Departement de Recherche Spatiale (DESPA) 
of the Observatoire de Paris and the Laboratoire d'Astrophysique of the Observatoire de Grenoble 
for supporting this work. I am especially grateful to Kent Wallace and the Spatial Interferometry 
Group at JPL, Jean-Luc Beuzit, Jerome Bouvier and the referee for giving comments. Data are 
based on observations obtained at the European Southern Observatory, La Silla, Chile 
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